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Abstract. We present here the results of a deep (130 ks) XMM-Newton observation of the cluster of galaxies 2A 0335-1-096. 
The deep exposure allows us to study in detail its temperature structure and its elemental abundances. We fit three different 
thermal models and find that the multi-temperature wdem model fits our data best. We find that the abundance structure of the 
cluster is consistent with a scenario where the relative number of Type la supernovae contributing to the enrichment of the 
intra-cluster medium is ~ 25%, while the relative number of core collapse supemovae is ~ 75%. Comparison of the observed 
abundances to the supernova yields does not allow us to put any constrains on the contribution of Pop III stars to the enrichment 
of the ICM. Radial abundance profiles show a strong central peak of both Type la and core collapse supernova products. Both 
the temperature and iron abundance maps show an asymmetry in the direction of the elongated morphology of the surface 
brightness. In particular the temperature map shows a sharp change over a brightness edge on the southern side of the core, 
which was identified as a cold front in the Chandra data. This suggests that the cluster is in the process of a merger with a 
subcluster. Moreover, we find that the blobs or filaments discovered in the core of the cluster by Chandra are, contrary to the 
previous results, colder than the ambient gas and they appear to be in pressure equilibrium with their environment. 

Key words. Galaxies: clusters: general - Galaxies: clusters: individual: 2A 0335-1-096 - cooling flows - intergalactic medium 
- Galaxies: abundances - X-rays: galaxies: clusters 


1. Introduction 

Clusters of galaxies are the largest known gravitationally bound 
structures in the universe. According to the standard cosmolog¬ 
ical scenario they form and grow along the filaments through 
merging with groups and individual galaxies. Optical and X- 
ray studies reveal that clusters of galaxies are still forming at 
the present epoch. 

The large effective area and superb spectral resolution 
of XMM-Newton together with the high spatial resolution of 
Chandra allow us to study clusters of galaxies with unprece¬ 
dented detail. Analysis of data obtained by these satellites led 
to a number of important results in the recent years. In the cen¬ 
tral parts of many clusters of galaxies the gas density is high 
enough that the radiative cooling time of the gas is shorter than 
the age of the cluster. As the gas cools the pressure decreases, 
which causes a net inflow toward the center of the cluster. Many 
clusters of galaxies indeed show a temperature drop by a factor 
of three or more with in the central 100 kpc radius (for a review 
on cooling flows see lFahianlll994 . Flowever, the spectra ob¬ 
tained by the XMM-Newton Reflection Grating Spectrometer 
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(RGS) show no evidence for strong cooling rates of gas be¬ 
low 30 - 50% of the maximum temperature of the ambient gas, 
which forces us to l ook for additional heating mechanisms in 
the cores of clusters llPeterson et al.l200ltlTamura et al.l2001bl: 
iKaastra et alj200ll) . The high resolution images from Chandra 
led recently to the discovery of cold fronts, associated with mo¬ 
tion of the cluster cores and to the identificat ion of fllamentary 
struct ure in the cores of a number of clusters ilMarkevitch et alJ 
1200^ . Spatially resolved spectroscopy of many clusters of 
galaxies sho ws a strongly central ly peaked distribution of metal 
abundances jTamura et al 12004 . 

Because of their large potential wells, clusters of galax¬ 
ies retain all the enriched material produced in the member 
galaxies. This makes them a unique environment for elemental- 
abundance measurements and for the study of the chemical en¬ 
richment history of the universe. 

In this paper we study the X-ray bright cluster 
2A 03354-096 using spatially-resolved and high-resolution 
spectra obtained during a 130 ks observati on with the 

S e an Photon Ir naging Camera (EPIC, iTurner et alJ 
IStriider et alJ 120011) and the Refl ection Grating 
Spectrometer (RGS, Iden Herder eTall 1200ih aboard XMM- 
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Newton jjansen et Ii] l200ll) . The properties of 2A 0335+096 
allow us to address here some of the above mentioned issues. 


2A 0335+096 was first detected as an X-ray source by 
Ariel V jCooke et alJll97^ . and was fo und to be associated 


with a medium comp act Zwicky cluster 


Schwartz et al 


np; 

1980ll. The p r esence of a cooling flow was 


hrst noted by iSchwartz et al ] ( Il980l) in the data oh ained by 
HEAP 1. X-ray obs ervations with EXO SAT dSingh et al J 19861 
Il988l) and Einstein dWhite et alJll99lh confirme d the presence 
of the cooling flow. Observa tions with ROSAT dSarazin et alJ 
Il992t lirwin & Sarazinl[i995l) show a hlamentary structure in 
the ce ntral region of the cooling flow. Observations with 
ASCA dKikuchi et al.lll99^ show a hint of a centrally peaked 
metallicity distribution . Usin g data obtained by BeppoSAX, 


De Grandi & Molendil d2Q0l[ l2002l) analyzed the metallic¬ 


ity and temperature profile of the cluster and found a cen- 
trally peaked me tallicity gradient. Using the same dataset, 
iF.ttori et ^d2002h estimated the total mass of the cluster within 
the region with overdensity of 2500 times the critical density to 
be ~ 1.6 X 10'"^ Mq, while the mass of the gas was found to be 
~ 2.0 X 10'3 Mq. 


Recent Chandra observation shows a complex structure 
in the core of the cluster: a cold front south of the center, 
unusual X-ray morphology consisting of a number of X-ray 
blobs and/or hlaments o n scales > 3 kpc, a long with two 
prominent X-ray cavities dMazzotta et al.l2003i) . Moreover, the 
Chandra observation shows that the cluster has a cool dense 
core and its radial temperature gradient varies with position 
angle. The radial metallicit y profile has a pronou nced central 
drop and an off-center peak dMazzotta et alj2003l) . A previous 
shorter observation with XMM-Newton shows an increase of 
the Fe abundance to ward the center with a strong central peak 
dTamura et alJl20()3) . The central galaxy of 2A 0335+096 is 
a cD galaxy with a very extended optical emission line region 
(Ha+[NII]) to the northeast of the galaxy. Moreover, the central 
region of t he galaxy is anomalously blue , indica t ing rec ent star 
formation dRomanishin & Hintzenll9^ . lEdgj d200lll reports 
a detection of CO emission (implying 2x10^ M© of molecular 
gas) and IRAS 60 pm continuum. These obseryations indicate 
a mass deposition rate of a cooling flow of < 5 M© yr“'. A 
radio study of 2A 0335+096 shows a radio source coincident 
with the central ga laxy, which is surrounded by a mini-halo 
dSarazin et alJl995l) . 

Throughout the paper we use Hq - 70 km s ' Mpc ', Q.m - 
0.3, Qa = 0.7, which imply a linear scale of 42 kpc arcmin”' at 
the cluster redshift of z = 0.0349. Unless specified otherwise, 
all errors are at the Icr confidence leyel. 


2. Observations and data reduction 

2A 0335+096 was obseryed with XMM-Newton on August 4th 
and 5th 2003 with a total exposure of 130 ks. The EPIC MOS 
and pn instruments were operated in Full Frame mode using 
the thin Alter. The exposure times of both MOS cameras and 
RGS were 130 ks and the exposure time of pn was 93 ks. 


2.1. EPIC analysis 

The raw data are processed with the 6.0.0 yersion of the XMM- 
Newton Science Analysis System (SAS), using the EPPROC 
and EMPROC tasks. We apply the standard Altering, for 
EPIC/MOS keeping only single, double, triple and quadruple 
pixel eyents (PATTERNS 12), while for EPIC/pn we make use 
of single and double eyents (PATTERNS 4). In both cases only 
eyents with FLAG==0 are considered. The redistribution and 
ancillary response files are created with the SAS tasks rmfgen 
and arfgen for each camera and region we analyze. 

In the spectral analysis we remoye all the bright X-ray 
point-sources with a flux higher than 4.8 x 10“*"' erg s“' cm“^. 

At the low energies, the EPIC cameras are currently not 
well calibrated below 0.3 keV. Eor pn we take a conseryatiye 
lower limit of 0.5 keV, below which, at the time of our analy¬ 
sis, there are still some uncertainties in calibration. At energies 
higher than 10 keV our spectra lack sufficient flux. Therefore, 
in most of the analyzed regions, our analysis of the MOS and pn 
spectra is res tricted to the 0.3 - 10 ke V and 0.5 - 10 keV range 
respectiyely. iRonamente et al 1 il200.5h showed that the system¬ 
atic uncertainties of the EPIC detectors in the energy range of 
0.3 to 1 keV are less than 10%. We include systematic errors to 
account for the uncertainties in the calibration and background 
in the spectral fit. T he applied systematics are adopted from 
iKaastra et alJ il2004ll . After a preliminary analysis we And that 
the pn instrument has a small gain problem. We minimize this 
effect by shifting the pn energy grid by 10 eV. The spectra ob¬ 
tained by MOSl, MOS2 and pn are then fitted simultaneously 
with the same model, while their relatiye normalizations are 
left as free parameters. 

To test whether our abundance measurements can be influ¬ 
enced by narrow band regions where the calibration of EPIC 
might still be problematic, we fit the spectrum of PKS 2155- 
304, a bright BE Lac object obseryed with thin Alter on 2004 
Noyember 23. The X-ray spectra of BL Lac objects are gen¬ 
erally well fitted by simple absorbed power-law, or broken 
power-law models. We And that the fit residuals at the posi¬ 
tion of the measured spectral lines are always consistent within 
the Icr uncertainty with the best fit absorbed power-law model. 
Thus, we can rule out the possibility that the measured abun¬ 
dances are seriously influenced by narrow band calibration un¬ 
certainties. 

2.1.1. Background modeling 

Since we want to inyestigate the cluster properties up to at 
least 9' from its core, we need a good estimate for the back¬ 
ground, especially at the dim outer parts of the cluster. The 
source Alls the entire held of yiew, which makes a direct mea¬ 
surement of the local background yery diffi cult. The commonly 
use d combined background eyent fists of ll .iimh et al .I (120021) 
and iRead & PonmanI (I2003]) . can be used in areas where the 
surface brightness of the source is high, or when the cluster 
has similar backgrou nd properties as the c ombined background 
fields. As sfiown by Ide Plaa et al.l (l2006l) relatiyely small sys¬ 
tematic errors in the background normalization can haye a big 
influence on the spectral fits in the outer parts of clusters. 
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Table 1. The CXB components used in the fitting of 
2A 0335+096. The power-law index of the extragalactic com¬ 
ponent is frozen to y = 1.41. The fluxes are determined in the 
0.3 - 10 keV band and, apart of the Local hot bubble, corrected 
for Galactic absorption using a column of Ah=2.5 x 10^' cm“^ 
(the local hot bubble is inside of the Galactic absorbing col¬ 
umn). 


Component 

kr (keV) 

Flux (erg s ‘ cm ^ deg ^) 

Local bubble 

0.082 

3.58 x 10-*^ 

Soft distant component 

0.068 

6.57 x 10-*'* 

Hard distant component 

0.127 

6.19 x 10^*^ 

Extra galactic power-law 


2.16 x 10^" 


The background can be divided into three main compo¬ 
nents: instrumental background, low energy particles (with en¬ 
ergy of a few tens of keV) accelerated in the Earth’s mag¬ 
netosphere (so called soft-protons) and the Cosmic X-ray 
Background (CXB). 

To minimize the effect of soft protons in our spectral analy¬ 
sis, we cut out time intervals where the total 10-12 keV count 
rate deviates from the mean by more than 3cr. The cleaned 
MOSl, MOS2 and pn event files have useful exposure times 
of 107 ks, 108 ks and 77 ks, respectively. 

To deal with the instrumental backg round we adopt a 
method developed bv Ide Plaa et alJ ll2006l) . As a template for 
the instrumental background we use data from closed filter 
observations, which we scale to the source observation using 
events detected outside of the field of view (the values of the 
scaling factors are 1.00 ± 0.02, 1.04 + 0.02 and 1.10 + 0.07 for 
MOSl, MOS2 and pn respectively). The scaling is performed 
by adding or subtracting a powerlaw with a photon index of 
y = 0.15 for MOS an d y = 0.24 for pn. F or detailed descrip¬ 
tion of the method see Ide Plaa et alJ ( l2006l) . 

We correct for the Cosmic X-ray Background (CXB) dur¬ 
ing the fitting. Since all the CXB photons enter the mirrors to¬ 
gether with the source photons, we can fit the CXB simulta¬ 
neously with the source spectra. In our models, w e use CXB 
components described bv iKuntz & Snowdenlll2000l) . They dis¬ 
tinguish 4 different components: the extragalactic power law 
(EPL), the local hot bubble (LHB), the soft distant com¬ 
ponent (SDC) and the hard distant component (HDC). The 
EPL component is made from the integrated emission of faint 
discrete sources, mainly from distant Active Galactic Nuclei 
(A GNs). According to one of the most recent determinations 
by De Luca & Molendil il2004l) the power law index of the EPL 
is 1.41 + 0.06 and its 2-10 keV flux is 2.24 ± 0.16 x 10^“ 
erg cm“^ s“' deg“^ (with 90% confidence). Since in our spec¬ 
tral analysis we cut out all the point sources with a flux higher 
than 4.8 x 10~^ "^ erg s~^ cm~^ we r educe the EPL flux in our 
extraction area. iMoretti et alJ ( l2003t) have made a compilation 
of number counts of X-ray point sources in two energy bands 
(0.5 - 2 and 2-10 keV) from a large source sample and deter¬ 
mined analytic formula for the number of point sources N(S ) 
with a flux higher then S. To determine, the contribution of the 
point sources brighter then our cut-off to the total flux of the 


EPL we calculate the integral S dS. We find, that 

the point sources we cut out make up 20 % of the total EPL 
flux. The 0.3 - 10 keV flux of the EPL in our model is thus 
2.16 X 10 ** erg cm“^ s * deg“^. The LHB is a local super¬ 
nova remnant, in which our Solar System resides. It produces 
virtually unabsorbed emission at a temperature of ~ 10® K. 
The SDC and HDC originate at a larger distance, they might be 
identified with the Galactic halo. Galactic corona or the Local 
group emission and are absorbed by almost the full Galactic 
column density. They have a temperature of 1 - 2 million K. 
We model each of the three soft background components (the 
LHB, SDC and the HDC) by a MEKAL model with tempera- 
tu res of 0.082 keV, 0.068 k eV, 0.127 keV respectively (based 
on iKuntz & Snowdeni2000l) . We determine the normalizations 
of these soft components by fitting the spectrum extracted from 
an annulus, with inner and outer radii of 9' and 12' respectively, 
centered at the core of the cluster. The contribution to the to¬ 
tal flux from cluster emission at these radii is ~ 60%, the rest 
of the emission comes from the CXB. To account for the clus¬ 
ter emission we use an additional thermal model. In our final 
spectral fits we fix the normalizations of the background com¬ 
ponents. The temperatures and the 0.3 - 10 keV fluxes of the 
background components are given in Table[2 

2.2. RGS analysis 

We extract the RGS spectra with SAS version 6.1.0 fo llowing 
the same method as described in iTamiira et in( l2001al) . In or¬ 
der to get spatial information we select the events from several 
rectangular areas on the CCD strip in the cross-dispersion di¬ 
rection. Because the cluster Alls the entire fleld-of-view of the 
RGS, we need a blank field observation to extract the back¬ 
ground spectrum. For this observation we choose a Lockman 
Hole observation with an effective exposure time of 100 ks. 
The flare subtraction is analogous to the method used with 
EPIC, but now we use the events from CCD 9 outside the cen¬ 
tral area with a cross-dispersion of \xdsp\ > 30" to make the 
lightcurve. This method was applied to both source and back¬ 
ground datasets. 

Because the RGS gratings operate without a slit, the result¬ 
ing spectrum of an extended source is the sum of all spectra in 
the (in o ur case) 5' x ~ 12' field of view, convolved with the 
PSF (see iDavislEooIl for a complete discussion about grat¬ 
ing responses). Extended line-emission appears to be broad¬ 
ened depending on the spatial extent of the source along the 
dispersion direction. In order to describe the data properly, the 
spectral fits need to account for this effect. In practice, this is 
accomplished by convolving the spectral models with the sur¬ 
face brightness profile of the source along the dispersion direc¬ 
tion jTamiira et alJl20o3) . For that purpose we derive for each 
extraction region the cluster intensity profile from MOS 1 in the 
0.8 - 1.4 keV band along the dispersion direction of RGS and 
we convolve this MOSl profile with the RGS response in or¬ 
der to produce a predicted line spread function (Isf). Because 
the radial profile for an ion can be different from the mean 
MOS 1 profile, this method is not ideal. Therefore we introduce 
two scale parameters, namely the width and centroid of the Isf. 
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These parameters are left free during spectral fitting in order 
to match the observed profiles of the main emission lines. The 
scale parameter s for the width is the ratio of the observed Isf 
width to the nominal MOS1 based Isf width. 


3. Spectral models 


For the spectral analysis we use the SPEX package 
( iKaastra et alj|l996b . We model the Galactic absorption us¬ 
ing the hot model of that package, which calculates the trans¬ 
mission of a plasma in collisional ionisation equilibrium with 
cosmic abundances. We mimic the transmission of a neutral 
plasma by putting its temperature to 0.5 eV. To find the best 
description of the cluster emission we fit several combina¬ 
tions of collisionally ionised equilibrium (CIE) plasma models 
(MEKAL) to the spectra: a single-temperature thermal model; 
a combination of two thermal models; a differential emis¬ 
sion measure (DEM) model with a cut-off power-law distri¬ 
bution of emission measures versus temperature {wdeni). The 
wdem model appears to be a good empirical approximation 
for the spectrum in cooling cores of c lusters of galaxies (e.g. 
iKaastra et alJ2Q0'^ de Plaa et alj2004l) . The emission measure 
Y - j rigHiidV (where and are the electron and proton 
densities, V is the volume of t he source) in the wde m model is 
shown in Eq. m adapted from iKaastra et alJ ll2004l) : 


^ _ f AT'/" 

dT ~\o 


Tmin <T < r^ax, 
elsewhere. 


( 1 ) 


The emission measure distribution has a cut-off at = 
cEniax- The cut-off c is set in this study to 0.1. Eor a oo 
we obtain a flat emission measure distribution. The emission 
measure weighted mean temperature Tmean is given by: 


r^max ry Jrp 

Jr„i„ dT ^ 

j ’fmax dY JT- 


( 2 ) 


By integrating this equation between Tmin and Tmax we obtain 
a direct relation between Tmean and Tmax as a function of a and 
c: 


_ (1 + l/a)(l 
“ (2+ l/a')(l -cF“+i) 


A comparison of the wdem m odel with the classic al cooling- 
flow model can be found in Ide Plaa et alJ ll2n0.5h . We note 
that the wdem model contains less cool gas than the classical 
cooli ng-flow model, which is co nsistent with recent observa¬ 
tions diPeterson et al.l20()ill2003l) . 

The spectral lines in the MEKAL model are fitted self con¬ 
sistently. Erom the fits we obtain the numbers of atoms of all 
elements with detected line emission. To see whether differ¬ 
ent elements show similar abundances with respect to solar, it 
is convenient to normalize these numbers with respect to the 
relative solar abundances. To make the comparison with pre- 
vious work easier, we use the solar abundances as given by 
lAnders & Greves^ ( ll989t) for this normali zation. The more re¬ 
cent solar abundance determinations (e.g. iGrevesse & Sauvad 
Il998t lLoddersll2003h give significantly lower abundances of 


O 



Fig. 1. The total spectrum of the core of the cluster (with a ra¬ 
dius of 3'), which contains the cooling core (upper panel) and 
the spectrum of the 3' -9' region (lower panel). The continuous 
line represents the fitted wdem model. The EPIC pn and MOS 1 
+ MOS2 spectra are indicated. 


oxyge n and neon than those measured by lAnders & Grevess'3 
dl989t) . Use of these new determinations would only affect the 
representation of the elemental abundances in our paper, but 
not the actual measured values, which can be reconstructed by 
multiplication of the given values with the normalizations. 

4. Global spectrum 

4.1. EPIC 

To compare the properties of the cluster in its cooling core re¬ 
gion with the properties outside of the cooling core region we 
extract two spectra. One is from a circular region with a radius 
of 3' centered on the X-ray maximum of the cluster and one 
is from an annulus with inner radius of 3' and outer radius of 
9'. The large temperature gradient and projection effects make 
it necessary to fit the spectra with a multi-temperature model. 
We use the wdem model, described in Sect. |3and fix the abun¬ 
dances of elements with weak lines, which we do not fit, to 
0.6 and 0.3 times solar in the inner and outer region respec- 
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Table 2. Fit results obtained by fitting the wdem model to the 
spectrum extracted from a circular region with a radius of 3' 
centered on the core of the cluster and from an annulus with 
inner and outer radii respectively 3' and 9'. For the Galactic 
absorption we use a value of Nii-2.5 x 10^' cm“^. Emission 
measures (F = J rignudV) are given in 10®® cm“^. Abundances 
are given with respect to solar. 


Parameter 

0-3' 

3-9' 

Y 

15.89 + 0.13 

6.53 + 0.01 

(keV) 

3.75 ± 0.02 

5.36 + 0.06 

a 

0.75 ± 0.02 

1.34 + 0.10 

kT 

mean 

2.64 ± 0.02 

3.46 + 0.05 

Mg 

0.57 ± 0.08 

0.21+0.06 

Si 

0.72 ± 0.03 

0.37 + 0.03 

S 

0.61 ± 0.03 

0.27 + 0.04 

Ar 

0.42 ± 0.05 

0.14 + 0.09 

Ca 

0.84 ± 0.06 

0.62 + 0.12 

Fe 

0.532 ± 0.007 

0.383 + 0.006 

Ni 

1.46 + 0.10 

0.77 + 0.15 

X-/d.o.f. 

948/525 

666/491 


Table 3. Abundance upper limits for elements with weak lines, 
which could not be reliably detected, determined from the 
fluxes at the expected line energies of their helium-like emis¬ 
sion. We also show the abundance of calcium calculated with 
the same method (it is consistent with the abundance deter¬ 
mined in MEKAL). 


Element 

I (phot m ^ s *) 

Abund. (solar) 

Ca 

0.252 + 0.019 

0.80 + 0.05 

Ti 

0.017 + 0.013 

1.5 + 1.1 

Cr 

0.020 + 0.010 

0.5 + 0.2 

Mn 

-0.003 + 0.009 

-0.2+ 0.5 

Co 

0.000 + 0.008 

0 + 2 



0.5 1 2 5 10 

Enerqy (keV) 


Fig. 2. Residuals of the fit to the EPIC total spectrum extracted 
from a circular region with a radius of 3', with the line emission 
put to zero in the model. 


tively. The high statistics of the spectra allow us to determine 
the temperature structure and elemental abundances of several 
elements in these regions very precisely. The best-ht parame¬ 
ters are shown in Table |2l 

We find that the core is cooler than the outer region. We also 
find that the outer part has a broader temperature distribution 
than the core. The spectra with the indicated important spec¬ 
tral lines are shown in Eig.[2 All the detected spectral lines are 
shown in Eig. |2 which shows the residuals of the EPIC spec¬ 
trum with line emission put to zero in the model, so all spec¬ 
tral lines are clearly visible. Eor the first time we see a feature 
at the expected energy of chromium in a cluster of galaxies, 
which corresponds to a detection with 2cr significance. Eirst, 
we attempt to measure the abundances of elements shown in 
Tabled 

Since oxygen is predominantly produced by type II super¬ 
novae, it is an important element for constraining the enrich¬ 
ment scenarios. Oxygen has strong H-like lines at 0.65 keV. 
The strong absorption toward 2A 0335+096 (see AnnendixiAl. 
uncertainties in the calibration of EPIC and relatively nearby 
iron lines do not allow us to constrain the oxygen abundance 
accurately by EPIC. However, the RGS with its high spectral 
resolution allows us to determine the oxygen abundance in the 
core of the cluster (see Tablel^. The 2p-ls neon lines at 1.02 
keV are in the middle of the iron L complex (lying between 
about 0.8 to 1.4 keV). The resolution of the EPIC cameras is 
not sufficient to resolve the individual lines in the iron L com¬ 
plex which makes the neon abundance determination by EPIC 
unreliable. However, the high resolution of RGS allows us to 
determine the neon abundance at least in the core of the cluster 
(see TablegJ. The K shell lines of magnesium at 1.47 ke V are 
also close to the iron L complex, which makes our magnesium 
abundance determinations by EPIC somewhat sensitive to our 
temperature and iron abundance model. The magnesium abun¬ 
dance determinations at larger radii, where the surface bright¬ 
ness of the cluster is relatively low might be influenced by the 
instrumental aluminum line at 1.48 keV. The silicon, sulfur, ar¬ 
gon and calcium lines lie in a relatively uncrowded part of the 
spectrum and their abundances are in general well determined. 
Iron has the strongest spectral lines in the X-ray band. At tem¬ 
peratures above 3 keV its Kcr lines are the strongest at about 
6.67 keV and 6.97 keV, while the iron L-shell complex ranging 
from about 0.8 keV to 1.4 keV dominates at lower tempera¬ 
tures. These spectral lines make the iron abundance determina¬ 
tions the most reliable of all elements. The nickel abundance is 
determined mainly from its K-shell line blends at 7.80 keV and 
8.21 keV, which are partially blended with iron lines. 

In Table |2l we see that the abundances in the outer part are 
always lower than in the core of the cluster. 

We also attempt to estimate the abundance upper limits 
for elements with weak lines (titanium, chromium, manganese, 
cobalt), which can not yet be fitted in MEKAL in a self- 
consistent way. We determine their abundances from the fluxes 
at the expected line energies of their helium-like emission (see 
Tabled. The feature at the expected energy of chromium cor¬ 
responds to a detection at a 2cr level, while the feature at the 
expected energy of titanium corresponds to a Icr detection. 
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Fig. 3. 1^' and 2"“^ order spectrum of 2A 0335+096 extracted 
from a 4' wide strip centered on the core. The continuous 
line represents the htted wdem model. Between 15 and 18 A 
Fe XVII and Fe XVIII emission lines are visible. On the x-axis 
we show the observed wavelength. 

Table 4. Fit results for the RGS spectra extracted from a 4' 
wide strip. The value for Vh is in 10^* cm“^. Iron is given with 
respect to solar and the other abundances with respect to iron. 


Parameter 

single-temp 

wdem-modA 

fVu 

3.14 + 0.03 

2.95 + 0.04 

kT (keV) 

1.80 + 0.02 


kT 

max 


3.36 + 0.13 

kT 

mean 


2.51+0.10 

a 


0.52 + 0.03 

Oabs 

0.47 + 0.02 

0.54 + 0.03 

N/Fe 

1.7+ 0.6 

1.3+ 0.4 

0/Fe 

0.55 + 0.05 

0.49 + 0.05 

Ne/Fe 

1.11 +0.10 

0.85 + 0.08 

Mg/Fe 

0.97 + 0.10 

0.97 + 0.08 

Fe 

0.524 + 0.018 

1.07 + 0.06 

Scale So 

0.59 + 0.12 

0.81 +0.14 

Scale Spe 

0.60 + 0.04 

0.73 + 0.04 

X^ / d.o.f. 

1552/790 

1173/788 


4.2. RGS 

In Fig.|3we show the hrst and second order RGS spectra ex¬ 
tracted from a 4' wide strip in the cross-dispersion direction of 
the instrument. The extraction region is centered on the core of 
the cluster. Despite the spatial broadening, the strong spectral 
lines of magnesium, neon, iron and oxygen are well resolved in 
both spectral orders. In Tablel^we list the ht results for these 
spectra, for both the single temperature and wdem model. 

The values of the hts show that the wdem model hts 
better than the single-temperature model. However, the x^ for 
the wdem is still high, probably because of the disturbed na¬ 
ture of the core of the cluster. There are some systematic dif¬ 
ferences between the two models. The abundances found ht- 
ting the spectrum with the wdem model are in general twice 
the single-temperature values and the line widths of oxygen 


and iron, indicated with the scale parameters sq and Sfe, are 
also higher for wdem. Because of these differences, we divide 
the abundances of the other elements, by the value for iron. 
Since the continuum is not well determined by the RGS, these 
ratios are more reliable than the absolute value. The relative 
abundances, given in units with respect to solar, for neon and 
magnesium are similar to iron, while the oxygen abundance is 
about half the value for iron. We are also able, using the excel¬ 
lent statistics, to derive the abundance for nitrogen. This value 
is also consistent with that for iron, although the error bars are 
large. The widths of the oxygen and iron lines, indicated by the 
scale values, are very similar in this cluster, which shows that 
the spatial distribution of oxygen and iron might be similar. 

Regarding the DEM temperature structure, the value of a 
(0.52 + 0.03) for wdem is smaller than the EPIC value of 0.75 + 
0.02 from within 3' from the core, probably, because of the dif¬ 
ferent extraction regions. The different extraction regions and 
the disturbed nature of the core (see sections and 0 make 
the abundances determined from the EPIC and RGS difficult to 
compare. The mean temperature kEmean of the wdem results is 
2.51 ± 0.10 keV, which is higher than the single-temperature 
ht. 

Be cause the solar oxygen abundance in lAnders & Grevessel 
(Il989ll is slightly overestimated, we need to free the oxygen- 
abundance in the Galactic absorption component (Gabs) in our 
spectral model. This ef fect in absorption was hrst observed by 
IWeisskopf et alJ ll2004l) in an observation of the Crab pulsar 
and later conhrmed in a cluster observation of Abell 478 by 
Ide Plaa et al 1 (l2004 . The values of 0.47 + 0.02 and 0.54 ± 0.03 
that we hnd for the oxygen-abundance in the Galactic absorp¬ 
tion com ponent Oah.. are slig htly lower than the RGS values 
found bv ide Plaa et alJ (120041) . but consistent with the expected 
value of 0.58 + 0.08 based on an updated value for the solar 
oxygen abundance bv lAllende Prieto et alJ(l200lh . 

5. Radial profiles 

We determine the radial temperature and abundance prohles of 
several elements using projected and deprojected spectra ex¬ 
tracted from circular annuli, centered on the X-ray maximum 
of the cluster. We use annuli with outer radii indicated on the 
top of Table |3 We also extract radial temperature and abun¬ 
dance prohles from RGS in the cross-dispersion direction. 

5.1. Projected spectra 

We ht the spectra extracted from the annuli with three different 
models, and compare the results. Due to high background in 
the 6th and 7th annulus we ignore the pn data at energies higher 
than 7.5 keV. In our model we hx the abundances of elements 
which we do not ht to 0.3 solar. The results are presented in 
Eig.l^and Tabled Eirst we ht the spectra with a single temper¬ 
ature model. We hnd that the model hts our data poorly and the 
X^ in the core of the cluster is unacceptably high. Based on the 
large x^ of the ht, in the central region we can discard the sin¬ 
gle temperature model. We therefore try to ht 2 thermal models 
simultaneously, with coupled abundances and with a hxed sep¬ 
aration between the two temperatures as T 2 = 2T\ (where Ti 
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Fig. 4. Comparison of fit results obtained by fitting the EPIC data with a single-temperature model, two-temperature model and 
the wdem model. 


and T 2 are the temperatures of the two thermal components). 
Fixing the temperature to ^2 = 27’i is a good first approxi¬ 
mation of the multi-temperature structure, since iKaastra et al.l 
ll2004 and IPeterson et al ]( |2003|) found that in almost all cases 
of their cluster samples, at each radius, there is negligible emis¬ 
sion from gas with a temperature less than one-third to half of 
the fitted upper temperature. Using the two temperature model 
our fits improve significantly. The improves most signifi¬ 
cantly in the core of the cluster, but we can see a slight im¬ 
provement also at outer radii. However, fitting the data with 
the wdem model further improves the x^ of our fits in 6 of 7 
extraction annuli. 


Fitting the data with a single-temperature thermal model we 
see that the abundances of all elements except argon drop in the 
core of the cluster and have a off-center peak. As we go from 
the single-temperature model to the two-temperature model 
and further to the wdem model, abundances of all elements in 
general increase, especially in the central regions of the cluster. 
The wdem model shows that most of the abundances peak in the 
core of the cluster, we detect a slight hint of a drop in the center 
only for calcium, silicon and magnesium. The dependence of 
the iron abundanc e on the ternperature model is known as th e 
Fe-bias (see, e.g. lRuotell2nn(i iMolendi & GastaldellolEnnl . 
As we go from one model to the other, the abundance value 
of magnesium varies the most. Due to its vicinity to the Fe-L 
complex its abundance is very sensitive to the temperature and 
iron abundance model. 


In general, all three models show a drop of abundances to¬ 
ward the outer parts of the cluster. 


5.2. Deprojected spectra 

To account for projection effects in the core of the cluster we 
extract spectra, deprojected under the assumption of a spheri¬ 
cal sym metry. We use the background event files of lUumb et alJ 
1120021) . We extract the deprojected spectra from 5 annuli in the 
inner 0' - 4', where the cluster is still bright enough and, as 
shown in Appendix the background subtraction does not 
have a significant influence on our results. The deprojected 
spectra represent the count rates from spherical shells cen¬ 
tered on the core of the cluster. Our extraction method and data 
analysis for the depr ojected spectra is described extensively in 
IKaastra et alJ ) l2004 . 

The fit results for the deprojected spectra are shown in 
Table The temperatures determined from deprojected and 
projected spectra are compared in Fig.|3 Since, in the process 
of deprojection we account for the hot plasma lying in front 
of the relatively cool core, the observed temperatures of the de¬ 
projected shells are lower than those determined without depro¬ 
jection. While the;^f^ of the fit of the single-temperature model 
improves significantly if we deproject the spectrum (compare 
with the x^ in Table 0, the single temperature thermal model 
still does not describe our spectra well. The wdem model de¬ 
scribes the deprojected spectra extracted from the core of the 
cluster significantly better than the single-temperature model. It 
shows, that the temperature gradients in the extraction regions 
and the intrinsic multi-temperature structure make it necessary 
to use multi-temperature models also when fitting deprojected 
spectra. Although the deprojected spectra give a better value 
for the temperature, with the deprojection we introduce more 
noise into the spectra so the determination of abundances be¬ 
comes more uncertain. 
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Table 5. Fit results obtained by fitting the EPIC data with a single-temperature model (1), two-temperature model (2) and the 
wdem model (3). Emission measures (F = J n^n-^dV) are given in 10^® cm“^. Abundances are given with respect to solar. 



Model 

0 - 0.5' 

p 

1 

in 

d 

1.0-2.0' 

2.0 - 3.0' 

d 

1 

d 

4.0 - 6.0' 

6.0 - 9.0' 

Y 

kT (keV) 
kT 2 (keV) 
(keV) 
(keV) 

a 

1 

2 

3 

1 

2 

3 

3 

3 

3.55 

3.42 

3.21 

1.87 + 0.01 

9 1 «4+0.005 
-0.002 
2.64 ± 0.02 

1.93 + 0.02 

0.58 + 0.01 

4.65 

4.82 

4.52 

2.32 + 0.01 
2.72 + 0.04 
3.14 + 0.02 

2.37 + 0.02 

0.48 + 0.02 

5.22 

5.45 

5.30 

2.78 + 0.01 

3.51+0.09 
3.95 + 0.03 
2.84 + 0.03 

0.65 + 0.02 

3.19 

3.31 

3.26 

3.04 + 0.01 
3.08 + 0.01 
4.39 + 0.04 

3.18 + 0.04 

0.62 + 0.04 

2.01 

2.08 

2.05 

3.18 + 0.02 

Q Q9+0.23 

4.81+0.05 

3.35 + 0.04 

0.79 + 0.04 

2.85 

2.75 

2.85 

3.14 + 0.02 

4.97 + 0.05 

3.34 + 0.08 

2.97 

2.56 

2.23 

3.32 + 0.04 
4.39 + 0.04 
5.84 + 0.14 

3.32 + 0.16 
8.02/‘3 

Mg 

1 

< 0.02 

< 0.18 

0.17 + 0.03 

0.07 + 0.05 

< 0.04 

< 0.01 

< 0.12 


2 

0.25 + 0.04 

0.14 + 0.05 

0.14 + 0.06 

< 0.13 

< 0.03 

< 0.01 

< 0.17 


3 

0.41 + 0.04 

0.43 + 0.06 

0.29 + 0.04 

0.26 + 0.06 

0.15 + 0.07 

< 0.07 

< 0.03 

Si 

1 

0.41 + 0.02 

0.58 + 0.03 

0.56 + 0.02 

0.46 + 0.03 

0.36 + 0.05 

0.12 + 0.04 

0.14 + 0.05 


2 

0.61 + 0.02 

0.62 + 0.03 

0.54 + 0.03 

0.42 + 0.04 

0.34 + 0.04 

0.14 + 0.04 

0.13 + 0.05 


3 

0.71 + 0.02 

0.76 + 0.03 

0.61 + 0.02 

0.53 + 0.03 

0.46 + 0.04 

0.23 + 0.03 

< 0.09 

S 

1 

0.32 + 0.04 

0.41 + 0.04 

0.38 + 0.02 

0.30 + 0.02 

0.20 + 0.06 

< 0.06 

< 0.03 


2 

0.58 + 0.03 

0.56 + 0.04 

0.40 + 0.04 

0.28 + 0.05 

0.22 + 0.06 

< 0.10 

< 0.04 


3 

0.67 + 0.03 

0.65 + 0.04 

0.51+0.03 

0.38 + 0.05 

0.35 + 0.06 

0.12 + 0.07 

< 0.10 

Ar 

1 

0.22 + 0.08 

0.16 + 0.08 

0.16 + 0.05 

< 0.07 

< 0.10 

< 0.05 

<0.03 


2 

0.52 + 0.06 

0.44 + 0.08 

0.21+0.10 

< 0.12 

< 0.20 

< 0.13 

<0.04 


3 

0.69 + 0.06 

0.50 + 0.09 

0.35 + 0.08 

0.19 + 0.12 

0.24 + 0.14 

< 0.18 

< 0.07 

Ca 

1 

0.51+0.12 

0.66 + 0.10 

0.85 + 0.07 

0.68 + 0.09 

0.53 + 0.15 

0.57 + 0.15 

0.27 + 0.21 


2 

0.66 + 0.09 

0.96 + 0.10 

0.92 + 0.12 

0.71+0.16 

0.66 + 0.14 

0.70 + 0.14 

0.30 + 0.23 


3 

0.74 + 0.09 

0.98 + 0.11 

1.02 + 0.10 

0.77 + 0.15 

0.74 + 0.17 

0.66 + 0.17 

0.70 + 0.23 

Fe 

1 

0.55 + 0.01 

0.61 + 0.01 

0.54 + 0.01 

0.45 + 0.01 

0.41 + 0.01 

0.32 + 0.01 

0.35 + 0.02 


2 

0.63 + 0.01 

0.57 + 0.01 

0.48 + 0.01 

0.40 + 0.01 

0.38 + 0.01 

0.29 + 0.01 

0.32 + 0.02 


3 

0.67 + 0.01 

0.62 + 0.01 

0.55 + 0.01 

0.43 + 0.01 

0.40 + 0.01 

0.31+0.01 

0.26 + 0.01 

Ni 

1 

0.68 + 0.11 

1.08 + 0.13 

0.57 + 0.08 

0.35 + 0.11 

0.37 + 0.20 

un. 

un. 


2 

0.89 + 0.09 

0.90 + 0.12 

0.53 + 0.13 

0.29 + 0.15 

0.22 + 0.17 

un. 

un. 


3 

1.65 + 0.10 

1.52 + 0.14 

0.99 + 0.10 

0.64 + 0.20 

0.69 + 0.17 

un. 

un. 

/ d.o.f. 

1 

2444/525 

1278/525 

1066/525 

774/525 

717/525 

909/499 

816/476 


2 

844/525 

753/525 

883/525 

680/525 

664/525 

789/499 

742/476 


3 

782/525 

664/525 

707/525 

596/525 

532/525 

607/499 

880/476 


-A- 


-tA- 


□ Deprojected - 
A Projected 


Radius (arcmin) 


Fig. 5. Comparison of the projected and deprojected mean tem¬ 
peratures obtained by fitting the wdem model. 


5.3. RGS radial profiles 

Because of the high statistics of the RGS spectrum shown in 
Fig.El we are also able to extract smaller strips in the cross¬ 
dispersion direction and make a radial profile up to 2' from the 
center. The results of these fits are shown in TableQ It confirms 
the radial temperature profile in the core found with EPIC, al¬ 
though the temperature determined with RGS is slightly higher. 
The reason for this discrepancy is probably the difference in 
extraction regions. The RGS spectra also contain a small con¬ 
tribution of photons originating outside the core radius due to 
the spatial extent of the source. Again, like in the full field- 
of-view case, the mean temperatures kT'mean are slightly higher 
than the single-temperature fits, although the values are consis¬ 
tent within the error-bars. The value for a is rather well con¬ 
strained in the region between -2' and 0.5' and shows a slight 
discontinuity around O', but the statistics do not allow to draw 
a conclusion. Prom the values it is clear that a DEM model 
{wdem) fits the RGS spectra better. 

Prom the RGS abundance profiles in the cross-dispersion 
direction, which are listed in Table 0 we confirm the spatial 
distribution of oxygen and iron indicated by the line widths 
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Table 6. Fit results of the deprojected EPIC spectra extracted from 5 annuli in the inner 4' and fitted with a single-temperature 
model (1) and wdem model (3). This region contains the cooling core, where the projection effects are the most important. 
Abundances are given with respect to solar. 


Parameter 

Model 

0 - 0.5' 

p 

1 

p 

1.0-2.0' 

2.0-3.0' 

3.0-4.0' 

Ah (10^‘cm-^) 

1 

2.36 ± 0.02 

2.37 + 0.02 

2.47 + 0.03 

2.48 + 0.05 

2.53 + 0.05 


3 

2.52 ± 0.06 

2.58 + 0.05 

2.64 + 0.08 

2.64 + 0.11 

2.66 + 0.16 

kT (keV) 

1 

1.65 + 0.02 

2.07 + 0.02 

2.68 + 0.02 

3.04 + 0.03 

3.11+0.05 

(keV) 

3 

2.36 ± 0.04 

2.72 + 0.04 

3.65 + 0.08 

4.28 + 0.13 

4.25 + 0.21 

kTmean (keV) 

3 

1.76 + 0.03 

2.09 + 0.04 

2.75 + 0.07 

3.10 + 0.12 

3.16 + 0.21 

a 

3 

0.52 + 0.03 

0.43 + 0.03 

0.49 + 0.04 

0.62 + 0.08 

0.53 + 0.14 


3 

0.6 + 0.2 

0.7 + 0.2 

0.3 + 0.2 

0.3 + 0.2 

< 0.5 

Mg 

1 

0.11+0.08 

0.34 + 0.10 

0.65 + 0.07 

0.35 + 0.13 

<0.16 


3 

0.55 + 0.10 

0.53 + 0.11 

0.76 + 0.14 

0.36 + 0.15 

<0.26 

Si 

1 

0.46 + 0.04 

0.73 + 0.05 

0.80 + 0.04 

0.73 + 0.07 

0.52 + 0.09 


3 

0.74 + 0.06 

0.84 + 0.06 

0.84 + 0.07 

0.71+0.08 

0.50 + 0.10 

S 

1 

0.42 + 0.04 

0.61+0.05 

0.52 + 0.04 

0.47 + 0.08 

0.19 + 0.12 


3 

0.66 + 0.06 

0.75 + 0.06 

0.61 + 0.07 

0.51+0.09 

0.25 + 0.10 

Ar 

1 

0.40 + 0.11 

0.50 + 0.11 

0.25 + 0.10 

< 0.14 

0.3 + 0.2 


3 

0.57 + 0.14 

0.69 + 0.12 

0.33 + 0.15 

< 0.21 

0.4 + 0.3 

Ca 

1 

0.64 + 0.20 

0.81+0.18 

1.00 + 0.13 

0.9 + 0.2 

0.4 + 0.3 


3 

0.41+0.20 

0.88 + 0.16 

1.06 + 0.18 

1.0+ 0.3 

0.5 + 0.4 

Fe 

1 

0.50 + 0.02 

0.65 + 0.02 

0.61 + 0.02 

0.46 + 0.02 

0.41 + 0.03 


3 

0.67 + 0.03 

0.66 + 0.03 

0.57 + 0.02 

0.45 + 0.02 

0.40 + 0.03 

Ni 

1 

0.62 + 0.17 

1.8+ 0.2 

1.53 + 0.18 

0.9 + 0.4 

0.6 + 0.5 


3 

1.3 +0.3 

2.0 + 0.3 

1.8+ 0.3 

0.8 + 0.4 

0.6 + 0.5 

/ d.o.f. 

1 

1090/587 

668/570 

892/587 

515/587 

520/587 


3 

622/587 

489/587 

526/587 

481/587 

515/587 


Table 7. Fit results for spatially resolved RGS spectra between 8-25 A. The fitted models are (1) single-temperature CIE and (3) 
wdem. We use 2cr upper limits. 


Parameter 

Model 

-2.0/-1.0' 

-1.0/-0.5' 

-0.5 / 0' 

0/0.5' 

0.5-1.0' 

1.0-2.0' 

kT (keV) 

1 

4.12 + 0.13 

3.37 + 0.10 

2.89 + 0.06 

2.53 + 0.05 

3.21 +0.12 

4.4 + 0.3 

kT 

max 

3 

6.6 + 0.3 

5.2 + 0.2 

4.6 + 0.2 

3.73 + 0.14 

5.3 + 0.8 

7.8 + 0.9 

kT 

mean 

3 

4.6 + 0.2 

3.70 + 0.15 

3.12 + 0.15 

2.72 + 0.11 

3.8+ 0.7 

5.1 +0.6 

a 

3 

0.79 + 0.08 

0.70 + 0.06 

0.94 + 0.10 

0.60 + 0.05 

0.7 + 0.4 

1.3 + 0.3 

O/Fe 

1 

0.28 + 0.07 

0.34 + 0.05 

0.40 + 0.04 

0.36 + 0.04 

0.35 + 0.07 

0.57 + 0.18 


3 

0.29 + 0.06 

0.34 + 0.05 

0.41 + 0.04 

0.36 + 0.04 

0.36 + 0.06 

0.47 + 0.14 

Ne/Fe 

1 

1.34 + 0.17 

1.34 + 0.15 

1.83 + 0.14 

1.54 + 0.12 

1.00 + 0.17 

1.4 + 0.5 


3 

0.60 + 0.15 

0.73 + 0.13 

0.98 + 0.10 

0.98 + 0.10 

0.56 + 0.15 

<0.3 

Mg/Fe 

1 

<0.4 

0.66 + 0.17 

0.40 + 0.12 

<0.4 

<0.5 

<0.2 


3 

0.5 + 0.2 

1.19 + 0.17 

1.09 + 0.12 

0.68 + 0.11 

0.42 + 0.18 

1.0 + 0.5 

Scale s 

1 

1.5 + 0.2 

0.98 + 0.10 

1.01 +0.10 

0.90 + 0.07 

0.77 + 0.13 

4.5 + 1.0 


3 

1.30 + 0.16 

0.78 + 0.10 

0.69 + 0.06 

0.62 + 0.05 

0.58 + 0.11 

3.0 + 0.7 

X- / d.o.f. 

1 

852/701 

977 / 701 

1491/701 

1337/701 

932/701 

923/701 


3 

757 / 698 

802 / 698 

800/698 

853 / 698 

851 /698 

839/698 


in the full RGS spectrum. The O/Fe ratio is consistent with a 
value between 0.3-0.4 over the whole profile. The neon and 
magnesium abundance tend to peak in the center slightly more 
than iron. Because the individual spectra in the RGS spatial 
profile have lower statistics than the total RGS spectrum, we 
use one scale parameter s for all measured lines. Therefore, the 
fitted scale parameter, which accounts for the broadening of the 
lines because of the spatial extend of the source, is the average 
of all line widths given in units of the width of cluster emission. 


These widths increase naturally toward the outer part of the 
cluster, because the spatial profile of the cluster also flattens. 

6. Temperature and iron abundance maps 

In order to investigate the spatial variations of the tempera¬ 
ture and metallicity we extract temperature and iron abundance 
maps of the 5.5'x5.5' region centered on the core of the cluster. 

First, we determine the temperature and iron abundance on 
a grid with bin-size of 30" x 30". For every bin we compute 
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Fig. 6. Temperature and iron abundance maps of the inner 5.5' x 5.5' region of the cluster. Panel a and b: temperature and iron 
abundance map with each 15" x 15" pixel htted with a single-temperature model. Panel c and d: mean temperature and iron 
abundance map with each 30" x 30" pixel fitted with a wdem model. Note that North is up. West is to the right. All maps have 
overplotted the same X-ray isophots. 


a redistribution and an ancillary hie. We ht the spectrum of 
each bin individually by the multi-temperature wdem model. 
The spectral hts are done with a hydrogen column density hxed 
to the global value 2.5 x 10^* cm“^ (see AnnendixlAl. The red- 
shift is hxed to the spectroscopically determined value of the 
central cD galaxy (z = 0.0349). The abundances of all ele¬ 
ments except iron in our model are hxed to 0.3 times the solar 
value (the exact value of these abundances does not inhuence 
the temperature and iron abundance determinations), while the 
iron abundance is left as a free parameter. We hx the param¬ 


eter a to 0.58, which is the value determined for the core of 
the cluster (see Table|5}. The free parameters are the maximum 
temperature Tmax, the iron abundance and the normalization of 
the wdem component. 

The temperature map, which shows the mean temperature 
Tmean of the wdem model, reveals an elongated temperature 
structure in the South-Southeast North-Northwest direction and 
a sharp temperature change over the brightness edge south of 
the core. The iron abundance map shows a highly centrally 
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Fig. 7. MOSl + MOS2 image of 2A 0335+096 with overplot¬ 
ted extraction regions. The concentric circles correspond to 
our extraction annuli used to investigate the radial profiles, the 
square shows the 5.5' x 5.5' region used to extract the tempera¬ 
ture and iron abundance maps and the small circles correspond 
to extraction regions 1, 2, 3 and 4 used to verify that the tem¬ 
perature and abundances vary with position angle. 


peaked iron abundance with an extension to the North (see the 
panel c and d in Fig.|^. 

High count rates in the 5.5' x 5.5' central region allow us to 
reduce the size of the bins to 15" x 15" (which is roughly the 
FWHM of the point spread function of the pn detector) and ex¬ 
tract a higher resolution temperature and iron abundance map 
fitting a single temperature model to each bin (see panel a and 
b in Fig. |^. The map shows a sharp temperature change on 
the southern side of the core and a large extension of the cool 
temperature gas to North-Northwest. In contrast to the wdem 
model fit, the iron abundance map extracted using a single¬ 
temperature model shows a central drop of iron abundance, 
with an abundance peak in a ’’high metallicity ring” around the 
X-ray core with an extension to the North. We verified that the 
central drop of iron abundance is evident also in a temperature 
map with 30"x30" bins fitted with a single temperature model. 
The central abundance drop is rather the result of the oversim¬ 
plified model for the temperature structure of the core of the 
cluster, than a real feature (see subsection l5.H . 

To verify the azimuthal differences in the temperature and 
the metallicity we extract spectra from 4 circular regions with a 
radius of 0.5' and fit them with the wdem model. The extraction 
regions are shown in Fig.0 The regions 1, 2 and 3 are centered 
at the same distance from the core (2') but at different position 
angles (North, West, South). Region 4 is at a distance of 2.75' 
from the core at South-Southeast. Apart from the wdem nor¬ 
malization, Tmax and the iron abundance we leave free and fit 
the parameter a. Ah and the silicon and sulfur abundance. The 
best-fit parameter values are given in Table |3 We confirm the 


Table 8. Parameter values obtained by fitting regions 1,2,3 and 
4 indicated in Fig. 0 with the wdem model. Emission measures 
(T = J n^n-ndV) are given in 10^® cm“^, temperatures are given 
in keV, Ah in 10^' cm“^ and abundances are given with respect 
to solar. 


Par. 

1(N) 

2(W) 

3(S) 

4 (SSE) 

Y 

4.55 

2.63 

1.84 

1.16 

Ah 

2.54 ± 0.05 

2.32 + 0.1 

2.49 + 0.09 

2.38 + 0.17 

kTniax 

3.18 + 0.10 

4.15 + 0.20 

4.2 + 0.25 

4.9 + 0.4 

^^mean 

2.69 + 0.13 

3.22 + 0.20 

3.25 + 0.24 

3.5+ 0.4 

a 

0.22 + 0.06 

0.41+0.10 

0.41+0.11 

0.6 + 0.3 

Si 

0.60 + 0.06 

0.70 + 0.10 

0.59 + 0.11 

0.50 + 0.24 

S 

0.46 + 0.08 

0.43 + 0.14 

0.17 + 0.14 

0.31+0.19 

Fe 

0.51+0.02 

0.39 + 0.03 

0.36 + 0.03 

0.37 + 0.04 

y / d.o.f. 

605/509 

512/509 

551/509 

598/509 


temperature and iron abundance asymmetry seen in the temper¬ 
ature maps. 

7. Properties of the core on smaller scales 

The high statistics of the data obtained during our deep XMM- 
Newton observation allow us to better study the thermal prop¬ 
erties of the blo bs/filaments found in th e core of 2A 0335+096 
with Chandra ilMazzotta et al.l I2OO3I) . We use a smoothed 
Chandra ACIS image in the energy band of 0.5 - 1.5 keV to 
locate the blobs/filaments. We create three extraction masks; 
one for the region of bright blobs/filaments-, one for the bright 
region around the blobs, but excluding the blobs (the ambi¬ 
ent gas)', and one for the bright core but excluding the com¬ 
plex structure (the envelope)', see Fig. |3 We use these masks 
to extract spectra from our EPIC-MOS data. The Pull Width 
Half Maximum (PWHM) of the on-axis Point Spread Punction 
(PSP) of MOSl and MOS2 at 1.5 keV is 4.3" and 4.4" respec¬ 
tively, while the PWHM of the on-axis PSP of EPIC-pn at the 
same energy is 12.5". Since the large PSP of EPIC-pn would 
cause a strong contamination of the spectra extracted from the 
region of blobs/filaments, we analyze only the data obtained by 
MOS. 

We fit the three spectra simultaneously. The envelope is fit¬ 
ted by a thermal model. The ambient inter-blob gas is fitted 
by two thermal models, while all the parameters of the sec¬ 
ond thermal model are coupled to the model with which we fit 
the envelope. The region of blobs/filaments is fitted with three 
thermal models, where the parameters of the second and third 
model are coupled respectively to the model we use to fit the 
ambient gas, and the envelope. This way we de facto deproject 
the temperature of the blobs. The only free parameters in the fit 
are the normalizations of the thermal components, the temper¬ 
atures, and the iron abundances. The values of the other abun¬ 
dances are fixed to the values determined by fitting the central 
region with wdem. The redshift is fixed to the redshift of the 
central cD galaxy. The best fit parameters are shown in Table|9] 
and the fitted spectra are shown in Fig.|^ The fit shows that the 
deprojected temperature of the blobs/filaments is 1.14 ± 0.01 
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Fig. 8. Smoothed Chandra image of the core of 2A 0335+096 
with overplotted areas used to extract spectra from the data ob¬ 
tained during our deep XMM-Newton observation. The spectra 
of blobs/filaments, ambient gas and the envelope were extracted 
respectively from the area indicated by black, white and gray 
contours. 

Table 9. Temperature and iron abundance of the 
blobs/filaments, the ambient gas around these features and of 
the envelope of gas surrounding the inner core. The superscript 
a indicates that the parameters are coupled. 



blo./fil. 

amb. gas 

envelope 

kT 

1.14 + 0.01 

1.86 + 0.01 

2.79 + 0.01 

Fe (solar) 

0.63“ ± 0.02 

0.63“ + 0.02 

0.61 + 0.01 


keV, while the temperature of the ambient gas is 1.86 + 0.01 
keV. The iron abundance of the blobs is coupled to the iron 
abundance of the ambient gas, while the iron abundance of the 
envelope is fitted separately and is found to be somewhat lower. 
Coupling the iron abundance of the blobs with that of the am¬ 
bient gas is necessary to avoid anti-correlation between the two 
parameters and it has no effect on the determined temperature 
difference between the two components. 

iMazzotta et alJ ( l2003l) assumed that the blobs are ellipsoids 
and selected eight r egions. We assume the regions and the lu¬ 
minosities given by IMazzotta et alJ ll2003l) and with our new 
temperature we calculate with a MEKAL model the density of 
the blobs. We find a value of ~ 8.5 x 10“^ cm“^. From our 
best fit emission measure of the thermal model of the ambi¬ 
ent gas we calculate the density of the gas and find a value of 
~ 5.5 X 10“^ cm~^. These density estimates and our new values 
of the temperatures of these two components therefore indicate, 
that the cold blobs/filaments and the ambient gas are in thermal 
pressure equilibrium (nkT ^0.1 keV cm“^). 



Fig. 9. Spectra of the blobs/filaments, the ambient gas around 
these features and of the envelope of gas surrounding the in¬ 
ner core. The vertical line helps to see that there is more Fe-L 
emission at low energies in the spectrum of the blobs. 



Fig. 10. The elemental abundances derived from a wdem fit of 
the inner 3' of 2A 0335+096, fitted by a linear combination of 
the yields of SN la, SN II and Pop III, for the SN la model 
W7 and Pop III model with a core mass of 130 Mq. The black 
empty circles indicate elements not used in the fit. 


8. Abundances and enrichment by supernova 
types la/ll and Population III stars 

We investigate the relative contribution of the number of su¬ 
pernova types la/II to the total enrichment of the intra-cluster 
medium (ICM). We also investigate the possibility of putting 
constraints on the contribution by Population III stars. We try 
to determine the relative numbers of different supernovae 


Nla.+ Nil + Npoplll 

where Asx is either the number of Type la supernovae (SNe la). 
Type II supernovae (SNe II) or Popupation III stars (Pop III) 
contributing to the enrichment of ICM. Since the progenitors of 
Type II supernovae and Type Ib and Ic (SNe Ib/Ic) are massive 
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Table 10. Relative numbers of SN la, SN II and Pop III contributing to the enrichment of the intra-cluster medium. We show 
the results of hts of linear combinations of the yields of SN la, SN II and Pop III stars, with two SN la yield models, to the 
abundances of 7 elements shown in Table|2determined for the core and the outer region of 2A 0335+096. For the core we show 
the results of fits with two different Pop III star models and the results of hts without Pop III stars. 


Model 

Type 

0' - 3' 3' - 9' 



Tlpopiii = 130 Mq ^Vd.o.f. 

Mpopiii = 65 Mq / d.o.f. 

no Pop III x^ / d.o.f. 

no Pop III x'^ / d.o.f. 

W7 

SN la 
SN II 
Pop III 

0.24 ± 0.04 

0.76 ± 0.04 36 / 4 

(4± 1) X 10^3 

0.18 + 0.01 

0.85 + 0.01 60/4 

(-2.7 + 0.4) X 10-3 

0.20 + 0.01 

0.80 + 0.01 70/5 

0.39 + 0.04 

0.61 + 0.04 44/5 

WDD2 

SN la 
SN II 
Pop III 

0.25 ± 0.04 

0.75 ± 0.04 35 / 4 

(-2 ± 1) X 10-3 

0.27 + 0.03 

0.72 + 0.03 43 / 4 

(6 + 18) X 10-3 

0.26 + 0.01 

0.74 + 0.01 43 / 5 

0.37 + 0.04 

0.63 + 0.04 8/5 



Atomic number 

Fig. 11. The elemental abundances derived from a wdem ht of 
the inner 3' of 2A 0335+096, htted by a linear combination of 
the yields of SN la and SN II. For SN la we use the WDD2 
model. The black empty circles indicate elements not used in 
the ht. 

stars, we associate the heavy-element yields from SNe II with 
those from SNe Ib/Ic. For nucleosynthesis products of SN II 
we adopt an average yield of stars on a mass range from 10 M© 
to 50 Mq calculated by iTsuiimoto et alJ ll 19951) assuming a 
Salpeter initial mass function. For n ucleosynthesis produc ts of 
SN la we adopt values calculated bv llwamoto et and 

investigate two different models (see Table llOt . The W7 model 
is calculated using a slow dehagration model, while the WDD2 
model is calculated using a delayed-detonation model. Note 
that delayed-detonation models are curTently favored over de¬ 
hagration models by the supernova comm unity, whereas the 
WDD2 is favored by llwamoto'et Ii]) ll999l) . The intergalactic 
medium might also have been signihcantly enriched by Pop III 
stars which exploded as pair instability supernovae. We inves¬ 
tigate two models for heavy e l ement yields of Pop III stars cal¬ 
culated bv iHeger & Woosl^ ll2002l) . with a core mass of 130 
Mq and 65 M©. 

We assume that the observed total number of atoms A, of 
the element i (determined from its abundance and the emission 
measure distribution) is a linear combination of the number 


of atoms Yi produced per individual supernova Type la (Y,ja), 
Type II (Y, u) and Population III star (Y,jii): 

Ni - aYjja + bYiii + cY/ju, (5) 

where a, b and c are multiplication factors representing the to¬ 
tal number of each supernovae that went off in the cluster and 
enriched the ICM. 

First we ht the abundances obtained for the cooling-core 
region (the inner 3') and determine the best-ht values for a, b 
and c, which we use to estimate the relative numbers of su¬ 
pernovae. We ht the abundance values obtained for the inner 
3' region shown in Table |2 The best ht values of the relative 
contributions are shown in Table 1 101 and Fig.^J In the Fig. 1101 
we also show our best measurement of the chromium abun¬ 
dance and the measured oxygen and neon abundances deter¬ 
mined by the RGS and normalized to the iron abundance de¬ 
termined by EPIC (the absolute abundances are not well deter¬ 
mined by RGS). Since the low signihcance makes the value of 
the chromium abundance uncertain and the abundance values 
of the oxygen and neon may be affected by systematic uncer¬ 
tainties arising from using different extraction regions for the 
RGS and EPIC, we do not include them in the htting. We see 
that all models are consistent with a scenario where the relative 
number of Type la supernovae contributing to the enrichment 
of the intra-cluster medium is 20 - 30%, while the relative num¬ 
ber of Type II supernovae is 70 - 80%. 

The 130 Mq Pop III star model fits show, that their relative 
number contributing to the enrichment is about two orders of 
magnitude lower than the number of SN la/II. However, they 
produce according to the models about two orders of magni¬ 
tude more mass per star than SN la/II so their contribution to 
the enrichment of ICM could still be signihcant. However, the 
relative numbers determined for Pop III are model dependent 
and the high values of the;^^^s do not allow us to conhrm their 
contribution. The 65 M© Pop III star model fits give unphysical 
values for their relative numbers. In Pig.[nwe show our best- 
fit model of the observed abundances as a linear combination 
of only SN la and SN II. This model gives us the same rela¬ 
tive number of SN la as the models containing Pop III stars. 
Also, using the WDD2 SN la model, the does not change 
with respect to models containing Pop III stars. The fit results 
show that Pop III stars are not necessary to explain the observed 
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abundance patterns, although their contribution to the enrich¬ 
ment of the ICM cannot be excluded. We conclude that this 
analysis does not allow us to put any constrains on the enrich¬ 
ment by Pop III stars. 

In Fig. and we see that the measured abundance of 
calcium is significantly higher than that predicted by the mod¬ 
els. The abundance of nickel is well fitted only using the W7 
model of Type la supernovae with Pop III stars. Fitting the 
abundances as a linear combination of SN la and SN II, the 
W7 model predicts more and the WDD2 model predicts less 
nickel than observed. Mainly the discrepancies between the 
model and the measured abundances of calcium and nickel are 
responsible for the large ;^f^s. The WDD2 model always pre¬ 
dicts significantly more chromium than the upper limit derived 
from our data-point. The best measured abundances of oxy¬ 
gen and neon are also not well reconstructed by the models. 
This may be either due to problems in the supernova models 
or due to the fact that they were determined from a different 
extraction region with a different instrument. We note that the 
involved uncertainties an d complications when interpre ting the 
observed abundances (see lMatteiicci & Chiar)ninil200-5li do not 
allow us to differentiate between the two SN la models. 

Since our fits of the core show that we cannot put constrains 
on the contribution of Pop III stars we fit the “outskirts” (3'-9') 
with models containing only SN la and SN II. The fits show a 
higher contribution of SN la in the outer region than in the cool¬ 
ing core region, however the large values and the involved 
systematic uncertainties do not allow us to make a firm state¬ 
ment. 

9. Discussion 

We have used three different thermal models in our analysis of 
2A 0335+096 and find that the multi-temperature wdem model 
fits the data best. We note that a good description of the tem¬ 
perature structure is very important to obtain correct values for 
elemental abundances. Fitting the EPIC and RGS data we de¬ 
termine the radial temperature and abundance profiles for sev¬ 
eral elements. The elemental abundances have a peak in the 
core of the cluster and show a gradient across the whole clus¬ 
ter. The abundance structure of the cluster is consistent with a 
scenario where the relative number of SN la is 20 - 30% of the 
total numb er of supernovae enriching the ICM. Contrary to the 
findings of iBaumgartner et*^ ll2005l) . the Population III stars 
are not necessary to explain the observed abundance patterns. 
However, their contribution to the enrichment of the ICM can¬ 
not be excluded. Extracting temperature and iron abundance 
maps we find an asymmetry across the cluster. We find that 
the deprojected temperature of the blobs/filaments in the core 
of the cluster is lower than the temperature of the ambient gas 
and the blobs appear to be in thermal pressure equilibrium with 
the ambient gas. Below, we discuss the implications of these 
results. 

9.1. Intrinsic temperature structure 

Our results show a presence of plasma emitting at a range of 
different temperatures within every extraction bin. In the core 


of the cluster, this may be caused in part by projection effects 
and by the strong temperature gradient. The multi-temperature 
structure we find around the central bin may also be caused 
in part by the large scale temperature asymmetry of the clus¬ 
ter, so in one extraction annulus we extract the spectrum of 
both the colder and of the warmer gas from North and South of 
the cluster center respectively. However, the non-isothermality 
detected by fitting the deprojected spectra cannot be caused 
only by the width of our extraction annuli in combination with 
the temperature asymmetry. The temperature gradients across 
our annuli and the temperature differences between North and 
South are not large enough for that. The results can be inter¬ 
preted as intrinsic multi-temperature stru cture at each position , 
which is in agreement with the results of iKaastra et alJ ll2004l) . 
Eurthermore, these results show that fitting the spectrum of the 
intra-cluster gas with a single-temperature thermal model is in 
general not sufficient and introduces systematic errors espe¬ 
cially in the determination of elemental abundances. 

9.1.1. Spatial temperature distribution 

The temperature distribution has a strong gradient in the core. 
In the inner 3' the deprojected temperature drops by a fac¬ 
tor of ~ 1.8 toward the center. In the outer part (3' - 9') the 
radial temperature distribution is flat. We do not see a tem¬ 
perature drop toward the outskirts. We And an asymmetry in 
the temperature distribution and iron abundance distribution. 
The temperature maps show the asymmetry in temperature dis¬ 
tribution in the direction of the elongated surface brightness 
morphology, roughly in the South-Southeast North-Northwest 
direction. The temperature drops faster toward the core from 
the South than from the North. The temperature change is the 
strongest over the co ld front South of the core identified by 
iMazzotta et alJ (120031) in Chandra data. The cold fronts are in¬ 
terpreted as boundaries of a dense, co oler gas moving throug h 
the hotter, more rarefied ambient gas dM arkevitch g t_a 1 ■||2000|) . 
Erom the pressure jump across the cold front IMazzotta et alJ 
(l2003h conclude that the dense central gas core moves from 
North to South with a Mach number of M ^ 0.75 + 0.2. The 
elongated, comet-like shape of the cold central core shown 
in our temperature maps supports this picture. The elongated 
shape of the cold gas core may hint that we see a surviving core 
of a merged subcluster moving through and being stripped by 
a less dense surrounding gas. The possible merger scenario is 
discussed in subsection l9.3l 

9.1.2. X-ray blobs/filaments 

IMazzotta et alJ (l2003l) found that the deprojected temperature 
of the blobs is consistent with that of the less dense ambient 
gas, so these gas phases do not appear to be in thermal pressure 
equilibrium. Therefore they propose a possibility of a signifi¬ 
cant unseen non-thermal pressure component in the inter-blob 
gas, possibly arising from the activity of a central active galac¬ 
tic nucleus (AGN). They also discuss two models to explain 
the origin of the blobs; hydrodynamic instabilities caused by 
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the motion of the cool core and “bubbling” of the core caused 
by multiple outbursts of the central AGN. 

Ho wever, the temperature measurements of iMazzotta et alJ 
1 I 2 OO 3 I) have large uncertainties due to the low statistics of 
the Chandra data. The superior spectral quality of our XMM- 
Newton data allows us to measure the temperature of the blobs 
much more accurately. In contrast to the previous hndings we 
hnd that the deprojected temperature of the region where they 
identihed the blobs/hlaments is colder by a factor of 1.6 than 
the temperature of the ambient gas and the blobs appear to be 
in thermal pressure equilibrium with the ambient gas. In the 
light of this new result a non-thermal pressure component in the 
inter-blob gas is not necessary to explain the observed blobs. 

9.2. Abundance distribution 

The high statistics of the data allow us to determine the abun¬ 
dances of the most abundant metals. For the hrst time we put 
constraints on the abundance of chromium, and upper limits on 
abundances of titanium, manganese and cobalt. The abundance 
distribution has a peak in the core of the cluster. The radial 
prohles show an abundance gradient toward the core across the 
whole analyzed region of the cluster. The abundance gradient is 
not restricted to the cooling-flow region. We observe centrally 
peaked abundance distributions for both SN la products (iron, 
nickel) and elements contributed mainly by SN II (magnesium, 
silicon). Contrary to our results fro m 2A 0335+096, coo ling 
core clusters (e.g. see the sample of iTamura et alJ ll2004l) l in 
general show a centraly peaked distribution of SN la products 
with respect to SN II products. The different abundance distri¬ 
bution in 2A 0335+096 might point toward a different enrich¬ 
ment history. In Sersic 159-03 the widths of the RGS lines also 
indicate that t he iron abundance is much more centrally peaked 
than oxygen (Ide Plaa et al.ll200^ . Contrary to Sersic 159-03, 
the RGS data in 2A 0335+096 indicate that the spatial distri¬ 
butions of oxygen and iron might be very similar. However, we 
have to be careful with interpreting the RGS results, because 
the strong temperature drop in the center of this cluster might 
cause that the RGS spectrum is dominated by the line emis¬ 
sion from the cool core and the line profiles of iron and oxygen 
follow the temperature structure more than the abundance dis¬ 
tribution of the elements. 

We estimate that the total mass of iron within the radius of 
126 kpc is ~ 2.1 X 10^ Mq. If all iron originates from SN la, 
we need ~ 2.8 x 10® SN la explosions to explain the observed 
amount of iron. Assuming 1.1 x 10'® years and 50 galaxies 
contributing iron via galactic winds and gas stripping we get an 
average SN la rate of one explosion per 200 years per galaxy. 

We show that the determined abundance values are sen¬ 
sitive to the applied temperature model. Especially single¬ 
temperature thermal models can introduce a bias in the abun¬ 
dance determination, resulting in a de tection of off-center abun¬ 
dan ce peaks reported for ex ample bv ISanders & FabianI (l2002i) 
and ijohnstone et all 1 I 2 OO 2 I 1 . 

The abundance maps show an asymmetry in the iron abun¬ 
dance distribution, with a higher abundance to the North of the 


cluster core. The abundance asymmetry further highlights and 
supports the picture of great dynamical complexity. 

9.3. Possibie merger scenario 

We observe an asymmetry of the temperature and iron abun¬ 
dance distribution across the cluster, along the surface bright¬ 
ness elongation, which points toward the poss i bility of an on¬ 
going merger in 2A 0335+096. iMazzotta et alJ(|2flfl3l) note that 
the orientation of the cold front and the asymmetry of the tem¬ 
perature distribution suggest that the cool core is moving along 
a projected direction that goes from North to South. However, 
it is not known whether the motion is due to a merger, or due to 
gas sloshing induced by an o f f-axis merger or by a flyby of an¬ 
other cluster. IMazzotta et aT ] (l2003l) also reason that if the cool 
core is indeed a merging subcluster, then in projection the po¬ 
sition of the subcluster is close to the cluster center. However, 
the impact parameter cannot be measured hence the subclus¬ 
ter may be passing through the cluster at any distance along 
the line of sight. Our detection of an asymmetry in the abun¬ 
dance distribution across the cluster, with a higher metallicity 
North of the core with respect to the South, supports the merger 
scenario, in which the subcluster had a higher metallicity. The 
most important enrichment processes, like ram pressure strip¬ 
ping of cluster galaxies and galactic winds, are not expected 
to produce the observed asymmetry in the spatial distribution 
of the iron abundance. In the dense cores of clusters of galax¬ 
ies the internal kinematics of the ICM is expected to quickly 
smooth out any existing inhomogeneities in the abundance dis¬ 
tribution. 

In Fig. ini we show an archival Hubble Space Telescope 
(HST) image of the core of the cluster with overplotted EPIC- 
MOS X-ray isophots. The HST image shows the central cD 
galaxy and a companion galaxy to the North-West (the bright 
source in the lower left is a foreground star). We see, that the 
X-ray brightness peak is offset to the South from the central 
cD galaxy, which also shows that the core is not i n a re laxed 
state. This was already shown hv IMazzotta et al ] (I 2 OO 3 I) who 
found that the centroid positions of the best fit double /3-model 
to the cluster surface brightness do not coincide with the posi¬ 
tion of the central cD galaxy, but are offset to the South-East 
and North-West. It is interesting to note, that the alignment of 
the central cD galaxy and the alignment of the central galaxy 
with the projected companion nucleus are following the same 
direction as the p ossible merger axis inferred from the X-rays. 
IPubinskil ( Il998l) has shown that the long axis of a brightest 
cluster galaxy always aligns with the primordial filament along 
which the cluster collapsed and with the long axis of the cluster 
galaxy distribution. If 2A 0335+096 formed along a filament 
aligned as the long axis of the central cD galaxy, we might be 
seeing the traces of the latest phase of the formation of the clus¬ 
ter, in which a subcluster is falling in along the same primordial 
filament. 

Although there is growing evidence for the merger sce¬ 
nario, with the available data we still cannot rule out an alterna¬ 
tive scenario of gas sloshing. Radial velocity measurements of 
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Fig. 12. An image of the core of the cluster obtained by 
the Hubble Space Telescope with overplotted X-ray isophots. 
North is up. West is right. The X-ray isophots are off-set from 
the central cD galaxy. The gray contour indicates the location 
of the cold blobs/filaments detected in X-rays. 


the individual galaxies in 2A 0335+096 are necessary to pro¬ 
vide the decisive evidence. 

9.4. Relative enrichment by supernova types la/ll and 
Population III stars 

All our fits show that the contribution of SN la is 20 - 
30% of the total num ber of supernovae enriching the ICM. 
iTsuiimoto et al.l(ll995l) found that the ratio of the total numbers 
(integrated over time) of supernovae Type la to Type II that best 
reproduces the Galactic abundances is Nu/Nu - 0.15, while 
the ratio determined for the Magellanic Clouds is 0.2 - 0.3, 
which correspond to relative numbers of Nia/Nia+ii = 0.13 and 
0.17 - 0.23 in the Galaxy and in the Magellanic clouds re¬ 
spectively. Supernova types la/II produce heavy-elements on 
different time scales during the chemical evolution of galax¬ 
ies, with SN II causing the enrichment usually in the early 
phases of galactic evolution and SNe la on a much longer 
time-scale, in the later phases of galactic evolution. The cur¬ 
rent relative frequencies of SNe la and II also depend on types 
of galaxies. According to the Lick Observatory Supernova 
Search (LOSS) about 42% of the supernovae observed in galax¬ 
ies within a redshift of z = 0.03 are SN la, while SNe II 
are rare in e arly-type galaxies and SNe la occur in all types 
of galaxies llvan den Bergh et ^ l2005h . The relative contri¬ 
bution to the enrichment of ICM by SN la in the cluster of 
galaxies 2A 0335+096 is between the Galactic value and the 
value of the current relativ e frequencies determined by LOSS 
( Ivan den Bergh et alJ2005l) . The value ofNiJN^+n = 20-30% 
is consistent with a picture of an early enrichment by SNe II 
material from star-burst galaxies and later contribution by 
Type la supernovae material. 


However, the scenario pr oposed by iTamura et al.l (12001 ai) 
and iFinoguenov et al ] ( l200ll) . where an early and well mixed 
SN II enrichment is followed by SN la enrichment in the cen¬ 
tral galaxy, produces over the Hubble time a centrally peaked 
distribution of SN la products relative to SN II products. As 
discussed in l9.2l we observe a centrally peaked abundance dis¬ 
tribution of both SN la and SN II products. Moreover, contrary 
to the expectations, comparing the relative enrichment in the in¬ 
ner and outer regions of the cluster the Mg/Fe, Si/Fe and S/Fe 
ratios point toward a higher contribution of SN la in the outer 
parts compared to the cool core. Such a difference in the en¬ 
richment history between the core and the outer parts would 
strongly support the scenario according to which the cool core 
is the surviving ICM of a merged subcluster tsee l9.3> . However, 
the involved uncertainties in the determination of the relative 

enrichment are high. _ 

As discussed by iMatteucci & Chiappinil (l2005l) the inter¬ 
pretation of the observed abundance ratios is not easy. The ob¬ 
served abundance ratios strongly depend on the star-formation 
history in cluster galaxies. Part of the supernova products might 
still be locked in the stars or in the potential wells of cluster 
galaxies. Therefore, the relative number of supernovae con¬ 
tributing to the enrichment of the ICM, which we infer from 
the abundance ratios, might be somewhat different from the rel¬ 
ative number of supernovae which exploded in the cluster over 
its life-time. 

The results also show that it is not necessary to include 
the Pop III stars in our models to explain the observed abun¬ 
dance patterns. Our magnesium, silicon, sulfur, argon and iron 
abundances can be reconstructed with a linear combination 
of yields of SN la and II for all tested SN la yield models. 
In the cluster 2A 0335+0 96 we cannot confirm the results of 
iBaumgartner et '^( l2005b claiming, that especially in the hot¬ 
ter clusters. Pop III stars are necessary to explain the observed 
silicon, sulfu r, calcium and argon abund ances. However, to test 
the results of IBaumgartner et d 1 (l2005t) an analysis of a large 
cluster sample including hotter clusters will be necessary. 

All fitted models give a significantly lower calcium abun¬ 
dance than we observe. This points toward uncertainties in the 
models of SN II yields, e.g. different explosion energies in the 
models can result i n a difference of an order of magnitude in 
calcium production dWooslev & WeaveTll995h . 

10. Conclusions 

We have analyzed spatially resolved and high resolution spec¬ 
tra of the cluster of galaxies 2A 0335+096 obtained during a 
deep XMM-Newton observation. We found that; 

- We unambiguously detect multi-temperature structure in 
the core of 2A 0335+096: for the determination of elemen¬ 
tal abundances, a multi-temperature model is mandatory. 
The wdem model is a good description of the temperature 
structure of the gas in 2A 0335+096. 

- The blobs/filaments found in the core of 2A 0335+096 are 
significantly colder than the ambient gas and they appear to 
be in pressure equilibrium with their environment. 
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- We detect a strong central peak in the abundance distribu¬ 
tions of both SN la and SN II products. 

- The relative number of SN la contributing to the enrich¬ 
ment of the intra-cluster medium in the central 130 kpc of 
2A 0335+096 is ~ 25%, while the relative number of SN II 
is ~ 75%. Comparison of the observed abundances to the 
supernova yields does not allow us to put any constrains on 
the contribution of Pop III stars to the enrichment of the 
ICM. We observe significantly higher calcium abundance 
than predicted by supernova models. 

- The detected asymmetry in the temperature and iron abun¬ 
dance distribution further supports the merger scenario, in 
which the subcluster had a higher metallicity. 
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up to 4' from the core and in the outer parts, where the back¬ 
ground starts to play a more important role the differences are 
large. This is shown in Fig. IA.2I and Table lA.ll We see that 
using the new background modeling the iron, silicon and mag¬ 
nesium abundances between 4' - 9' from the core are lower 
than those derived using the standard background event hies. 
The large magnesium abundances at the outer radii obtained us¬ 
ing the standard background event hies are due to the incorrect 
subtraction of the instrumental aluminum lines. Abundances of 
other elements are consistent within the errorbars. 


Appendix A: The influence of the background 
subtraction on radiai profiies 

We compare the radial temperature and abundance pro- 
hles obtained using th e combined background event lists of 
Read & PonmanI ll200.3ll with those determined using the back¬ 


ground model described in subsection 12.1.11 We compare the 
values obtained using the wdem model, which was shown to 
describe the spectrum in the cool cores of clusters well (see 
Sect. 0. We derive a radial hydrogen column density prohle, 
maximum temperature prohle, a prohle and abundance prohles 
of several elements: magnesium, silicon, sulfur, argon, calcium, 
iron and nickel. 

Fitting the data u sing the background event hies of 
Read & Ponman ( 200.31) we hnd an average hydrogen column 
density of (2.49 + 0.03) x 10^' cm“^, which is signihcantly 
higher than the Galactic value of 1.80 x 10^* cm~^ deter¬ 
mined from the 21 cm radio emission dPickev & Fockm^ 
Il990l) . This high absorption column density was fou nd before 
in ROSAT PSPC spectra by llrwin & SarazinI ill99.5ll and in a 

[ irevio us. shorter XMM-Newton observation bv iKaastra et al.l 
2 OO 3 I) . The additional absorption might be caused by Galactic 
dust clouds or molecular clouds. The measured radial absorp¬ 
tion distribution has a slight peak between 2' and 3', but the 
errobars are large and it is still consistent with a hat distribu¬ 
tion (see Fig. lA.lli . Therefore, when htting the data with our 
background model, we hx the absorption to 2.5 x 10^' cm“^. 

The left top panel in Fig. 10 shows that the radial tem- 
perature prohles derive d with the background event hies of 
iRead & PonmanI ) l2003l) and with our background model are 
consistent with each other up to 6' from the core of the clus¬ 
ter. The best ht abundance values are, however, consistent only 
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Fig.A. 2. Comparison of fit rem its obtained by fitting the wdem model to EPIC spectra using the combined background event 
files of Read & PonmanI il2003li and using our new background subtraction method. 
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Table A.l. Fit results o btained by fitting the wdem model to EPIC spectra using (R) the combined background event lists of 
Read & PonmanI ( l2003l) and using (M) the new background subtraction method. Abundances are given with respect to solar. 


Parameter 

backg. 

o 

1 

o 

p 

1 

d 

1.0-2.0' 

2.0 - 3.0' 

b 

1 

p 

4.0 - 6.0' 

6.0 - 9.0' 

Nh (lO^'cm-^) 

R 

2.48 ± 0.03 

2.50 + 0.03 

2.52 + 0.03 

2.55 + 0.03 

2.51+0.04 

2.48 + 0.03 

2.40 + 0.05 

(keV) 

R 

2.66 ± 0.02 

3.16 + 0.02 

3.98 + 0.04 

4.40 + 0.06 

4.73 + 0.07 

5.04 + 0.10 

6.28 + 0.23 

(keV) 

M 

2.64 ± 0.02 

3.14 + 0.02 

3.95 + 0.03 

4.39 + 0.04 

4.81+0.05 

4.97 + 0.05 

5.84 + 0.14 

(keV) 

R 

1.95 ± 0.02 

2.38 + 0.02 

2.92 + 0.03 

3.17 + 0.06 

3.35 + 0.07 

3.49 + 0.09 

3.81 + 0.25 

fcTniean (keV) 

M 

1.93 ± 0.02 

2.37 + 0.02 

2.84 + 0.03 

3.18 + 0.04 

3.35 + 0.04 

3.34 + 0.07 

3.32 + 0.16 

a 

R 

0.57 ± 0.01 

0.49 + 0.02 

0.57 + 0.02 

0.64 + 0.05 

0.71+0.06 

0.83 + 0.07 

1 4+0.8 

a 

M 

0.58 ± 0.01 

0.48 + 0.02 

0.65 + 0.02 

0.62 + 0.04 

0.79 + 0.04 

1 Q 1+0.23 
^•^^-0.02 

2.37/f 


M 

0.95 ± 0.07 

0.79 + 0.10 

0.16 + 0.06 

0.1+0.08 

< 0.03 

< 0.01 

<0.01 

Mg 

R 

0.49 ± 0.06 

0.45 + 0.06 

0.45 + 0.07 

0.30 + 0.08 

0.22 + 0.08 

0.41 + 0.09 

0.57 + 0.11 


M 

0.41 ± 0.04 

0.43 + 0.06 

0.29 + 0.04 

0.26 + 0.06 

0.15 + 0.07 

< 0.07 

< 0.03 

Si 

R 

0.76 ± 0.03 

0.78 + 0.03 

0.68 + 0.03 

0.54 + 0.04 

0.48 + 0.04 

0.39 + 0.04 

0.37 + 0.05 


M 

0.71 ± 0.02 

0.76 + 0.03 

0.61 + 0.02 

0.53 + 0.03 

0.46 + 0.04 

0.23 + 0.03 

< 0.09 

S 

R 

0.70 ± 0.03 

0.68 + 0.03 

0.53 + 0.03 

0.40 + 0.04 

0.35 + 0.05 

0.20 + 0.05 

< 0.15 


M 

0.67 ± 0.03 

0.65 + 0.04 

0.51+0.03 

0.38 + 0.05 

0.35 + 0.06 

0.12 + 0.07 

< 0.10 

Ar 

R 

0.66 ± 0.06 

0.55 + 0.07 

0.35 + 0.08 

0.20 + 0.10 

0.32 + 0.12 

< 0.21 

< 0.04 


M 

0.69 ± 0.06 

0.50 + 0.09 

0.35 + 0.08 

0.19 + 0.12 

0.24 + 0.14 

< 0.18 

< 0.07 

Ca 

R 

0.74 + 0.10 

0.97 + 0.11 

0.98 + 0.13 

0.79 + 0.17 

0.78 + 0.18 

0.53 + 0.15 

< 0.25 


M 

0.74 ± 0.09 

0.98 + 0.11 

1.02 + 0.10 

0.77 + 0.15 

0.74 + 0.17 

0.66 + 0.17 

0.70 + 0.23 

Fe 

R 

0.67 ± 0.01 

0.63 + 0.01 

0.52 + 0.01 

0.44 + 0.01 

0.40 + 0.01 

0.36 + 0.01 

0.37 + 0.02 


M 

0.67 ± 0.01 

0.62 + 0.01 

0.55 + 0.01 

0.43 + 0.01 

0.40 + 0.01 

0.31+0.01 

0.26 + 0.01 

Ni 

R 

1.70 + 0.13 

1.60 + 0.14 

0.73 + 0.16 

0.73 + 0.18 

0.76 + 0.20 

0.94 + 0.22 

1.9 + 0.3 


M 

1.65 + 0.10 

1.52 + 0.14 

0.99 + 0.10 

0.64 + 0.20 

0.69 + 0.17 

un. 

un. 

A- / d.o.f. 

R 

726/525 

764/525 

814/525 

721/525 

653/525 

712/499 

1281/499 


M 

782/525 

664/525 

707/525 

596/525 

532/525 

607/499 

923/499 






















